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A? 

Abstract: A detailed study of the thermodynamics of the 0(N = 3) model in 1+1 
dimensions is presented, employing a 2PI resummation prescription as well as fully non- 
perturbative finite-temperature lattice simulations. The analytical results are computed 
using the CJT formalism and the auxiliary field method to one- and to two-loop order. 
The lattice results are obtained via the integral method using a Monte Carlo simulation 
for various lattice spacings. The analytical and lattice results for pressure, trace anomaly, 
and energy density, resembling closely those of four-dimensional Yang-Mills theories, are 
compared with each other. We find that to one-loop order there is a good correspondence 
between the CJT formalism and the lattice study for low temperatures. However, at high 
T the two-loop calculation fares better, correcting for the overestimation from the former 
approximation. 

Keywords: 0(3) model, CJT formalism, lattice thermodynamics, low-dimensional models 



Contents 



1 Introduction 1 

2 The model and the analytic results 3 

2.1 Generating functional 3 

2.2 The CJT effective potential 4 

2.3 One-loop approximation 5 

2.4 Regularisation 7 

2.5 Renormalisation 8 

2.6 Two-loop approximation 9 

3 Lattice simulation 12 

3.1 Thermodynamical observables 12 

3.2 Scale setting 13 

3.3 Numerical results 14 

3.3.1 Beta-function 14 

3.3.2 Thermodynamics 14 

4 Comparison of analytical and lattice results 16 

5 Summary 19 
A Algorithms and autocorrelations 20 



1 Introduction 

As discussed by Polyakov [1], the two-dimensional O(N) nonlinear sigma model has inter- 
esting features in common with four-dimensional non-Abelian gauge theories. For these 
reasons it has been investigated in the past as a suitable toy model for QCD, see Refs. [2-7] 
and refs. therein. 

In 1 + 1 dimensions, the coupling constant is dimensionless, the theory is renormalizable 
both perturbatively and in the 1/N expansion, and has also a negative beta- function, thus 
showing asymptotic freedom [1]. Moreover, the two-dimensional O(N) model is conformally 
invariant at the classical level, but, just as Yang-Mills theories in four dimensions, at the 
quantum level an energy scale emerges due to quantum corrections (trace anomaly). In 
fact, it exhibits a nonperturbative mass gap m = /xexp(— 2ir/ g 2 ) generated dynamically due 
to the interactions, where g is the coupling constant and fi the renormalisation parameter. 
Note that: (i) since the mass is non-analytic for g — > 0, it vanishes in perturbation theory; 
(ii) unlike the case of more than two dimensions, here no spontaneous breaking of the 
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global O(N) symmetry takes place: this is related to the Mermin-Wagner-Coleman theorem 
[8, 9], which forbids spontaneous breaking of a continuous global symmetry in a (1 + 1)- 
dimensional homogeneous system at any temperature. 

The effective potential at zero temperature has been investigated in Refs. [10, 11]. In 
Ref. [12] the model was studied at finite temperature in perturbation theory and to leading 
order in the large- N limit; this was subsequently extended in Refs. [6, 7] to next-to-leading 
order (NLO) in the 1/N expansion. It is found that the NLO effective potential contains 
temperature-dependent divergences. Fortunately, it turns out that the effective potential 
can still be renormalised at the minimum. 

In this work we study the nonlinear 0(N = 3) model in 1+1 dimensions at nonzero 
temperature both analytically and numerically. On the analytical side, we compute the 
thermodynamic properties by employing the CJT (Cornwall-Jackiw-Tomboulis) formalism 
[13] and using the auxiliary field method in the one- and two-loop approximation. We 
present analytical expressions for the renormalised pressure, the trace anomaly, the energy 
density, and the quasi-particle mass of the scalar fields. On the numerical side, we compute 
the thermodynamic quantities via "first-principles" Monte Carlo lattice calculation using 
the integral method. We then compare the analytic and the lattice results: at one loop, 
a good agreement is found for small temperatures, but the analytic result for the pressure 
is too large when T increases. When going to the two-loop approximation, a considerable 
improvement in the high-T domain is obtained. 

Furthermore, for N = 3 the model exhibits instantonic solutions both at zero [2] and 
nonzero temperature [14, 15]. In fact, the 0(3) model in 1 + 1 dimensions is topologically 
equivalent the CP 1 model, whereas CP n ~ l theories admit instantons for each value of n 
[16, 17]. The instanton solutions are topologically nontrivial field configurations localised 
in space and time, associated to classical, finite-action solutions of Euclidean field equa- 
tions. They give rise to tunneling processes between different classical vacua [18]. For this 
reason they have been regarded in the framework of Yang- Mills theories SU(N > 1) (see 
for instance [19, 20] and refs. therein) as a possible crucial ingredient for understanding 
nonperturbative properties, such as confinement and the emergence of a mass gap. 

Finally, one may wonder what is the role of instantons and calorons in the thermody- 
namics of this model; although our analytical method does not distinguish topologically 
trivial and nontrivial field configurations, we suspect, in agreement with [21], that the role 
of instantons and/or calorons can be relevant. 

The manuscript is organised as follows: in Sec. 2 we present the model, derive the 
effective potential to one- and two-loop order, and derive the analytic expressions for pres- 
sure, trace anomaly and energy density. In Sec. 3 we present a detailed lattice study of the 
model employing the integral method. In Sec. 4 we show and discuss the comparison of 
analytic and lattice results. The reader interested only in the results may skip the details 
of Sees. 2 and 3 and go directly to Sec. 4. Finally, conclusions and outlooks are described 
in Sec. 5. 

Throughout the paper, we will use natural units K= c = fee = 1. All computations at 
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nonzero temperature T are done in the imaginary-time formalism, with the shorthand 
where we have used the Euclidean momentum K = (k x = k, k T ). 



XI) 



2 The model and the analytic results 
2.1 Generating functional 

The nonlinear O(N) model in 1 + 1 dimensions at nonzero temperature is defined by the 
following generating functional: 



Z = M / £><£> 5 



N 

7/2 



exp 



l/T 



dT 



dx Cq 



(2.1) 



where g is the dimensionless coupling constant. The N fields of the model are contained 
in the column matrix <I>, which for future purposes we write as 



¥ = (o-,7ri,...,7rjv-i) 



(2.2) 



Moreover, we also use the notation $ 2 = and 7?* = (7Q, . . . , vrj\r_i) . 

The Lagrangian Cq entering Eq. (2.1) is a simple Lagrangian for N massless scalar 
fields, which in Euclidean space reads: 



(2.3) 



with d\ = d x , 82 = d T . 

Due to the 5-function entering the partition function (2.1) the fields of the model are 
constrained by the condition 



N 



(2.4) 



This nonlinear constraint enforces the N degrees of freedom of the model on an (N — 1)- 
dimensional hypersphere and induces interactions between them. Using the mathematically 
well-defined (i.e. convergent) form of the usual representation of the functional <5-function 
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J 2 



lim M \ Va exp < — / 

£^o+ J J 



l/T r+oo 
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J -oo 
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(2.5) 

the generating functional and the corresponding Lagrangian can be rewritten as follows: 



Z = lim M I Va V<5> exp 

0+ 



/oo 
dxC 
-oo 



(2.6) 
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where 



C 



a { $ 2 



(2.7) 



In the last expression a is an unphysical auxiliary field serving as a Lagrange multiplier. In 
the next Subsections we shall use the form in Eqs. (2.6) and (2.7) for our explicit calculation 
of the thermodynamical properties. Interestingly, by further eliminating a we obtain the 
following equivalent form of the generating functional: 

l<S> 2 -N/g 2 y 





rl/T /• oo 




Z = lim + N J V$> exp < 


— / dr / dx 






JO J-oo 



4e 

in which the "infinitely steep" mexican-hat form of the potential is evident. 
2.2 The CJT effective potential 

In order to compute thermodynamic quantities we use the CJT formalism [13] which is 
a self-consistent resummation prescription to compute the effective potential for a given 
theory. 

The first step consists in shifting the fields a and a. by their nonvanishing vacuum 
expectation values, 

a— >4> + cr , a — > Qo + a . (2-9) 
As a consequence, the Lagrangian (2.7) takes the form 

C = \d^¥d^ + ^(a + a) 2 + % -(a + a) {a 2 + vf 2 + 2<t0 + (f 2 - . (2.10) 

As one can see from Eq. (2.10) this shift produces a bilinear mixing term, iaa<f>, which 
renders the propagator matrix non-diagonal in the fields a and a. There are two possible 
ways to handle the ensuing mixing. One is to keep this term and allow for a non-diagonal 
propagator which mutually transforms the fields a and a into each other [22, 23]. The 
second way, used in this work and based on the study of Ref. [24], consists in performing 
the shift 

a — >a-4— a , (2.11) 

iVe 



which eliminates the mixing term. In Ref. [24] it was explicitly demonstrated up to two- 
loop order that these two methods to handle the mixing term are equivalent, i.e. they 
yield the same results for the effective potential and for the gap equations. However, the 
shift introduced in Eq. (2.11) considerably simplifies the calculations. From the resulting 
Lagrangian, 



r 1 aa , 1 a a , a ( ■ , <t> \ , n i I- \ , a 
L = -d^ad^a + -0^0^ + — Ua + — I + — {ta ) + — e 



+-a(a 2 + vr 2 ) + ^-(a 2 + vr 2 ) + ~a 



N 
9 



~2 J + 2 a ° ' 



(2.12) 
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one can immediately deduce the inverse tree- level propagators and the tree- level masses: 

Dr l (K;<t>,a Q ) = K 2 + mf , i = a, tt; (2.13) 

ml = ia + — 5 m l = i^o I ( 2 -14) 
D- 1 = ml = e . (2.15) 

Within the aforementioned CJT formalism the standard expression for the effective 
potential is given by 

Vett = V e fi{4>, ao; G a , G n , G a ) = 

U{4>,a Q )+ \ I [lnGr 1 (K) + Dr 1 (K;<p,a )G l (K)-l} + V 2 , (2.16) 

. - Jk 

where U(4>, ao) is the tree-level potential just as in Eq. (2.7), Di(K; (j), ao) are the tree-level 
propagators, Gi{K) are the full propagators in momentum space and V2 contains all two- 
particle-irreducible (2PI) self-interaction terms. The full propagators and the condensates 
are evaluated by using the stationary conditions for the effective potential: 

5V ef[ 5V eS 5V cS n . , n 1 ^ 

" ° ' 7a7 = ° ' 6G-(K) =0 ; l = ^ a ■ ( 7) 

For the full propagators one obtains 

Gr 1 {K) = K 2 + M 2 ■ M 2 =m? + S i , (2.18) 

where Sj denotes the self-energy 

E - = 2 mK) ■ < 2 - 19) 

2.3 One-loop approximation 

Restricting to one- loop order, the 2PI contribution to the effective potential is equal to 
zero, V2 = 0. Thus, the effective potential is given by: 

V cS = l -aoU 2 -^\+ e -a 2 +\ Y, j [^Gr\K)+Dr\K; ( f>,a )G i (K)-l) . (2.20) 

i=cr,TT,a 

Using the stationary conditions (2.17) one derives the following equations for the two 
condensates: 



() -<>[ inn + 4z [ K G ^ ! ; (2 ' 2I) 



Ne 



(2.22) 



b 2 -^ + I G a (K) + {N - 1) / G n (K) , 
9 Jk Jk 

and for the full propagators: 

Gr\K) = K 2 + M 2 =Dr\K-^,ao) ; (2.23) 
M 2 = mf + Si = mf . (2.24) 
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In fact, since V2 = 0, the self-energy Sj vanishes in the one-loop approximation. Substitut- 
ing the right hand side of Eq. (2.22) for icto gives the following equations for the physical 
condensate (ft and for the masses: 

= ^{ m * + mI k G(7{k) ) 5 (2 - 25) 

Ml =Ml + ^ ; (2.26) 



M 



2 



2 



,2 



[ G a (K) + (N-l) [ G n (K) 
9 Jk Jk 



(2.27) 



We now study separately the cases (ft 7^ and (ft = and show that, in the present two- 
dimensional model, the latter holds. 

(i) We assume (ft 7^ (corresponding to spontaneous symmetry breaking) . When taking 
the limit e — > + the degree of freedom denoted by a becomes frozen due to an infinitely 
heavy mass, see Eq. (2.26): 

lim Ml = lim ^t- = 00 , (2-28) 



from which it follows that 



- / G a (K) = . (2.29) 
e Jk 

Moreover, from Eq. (2.25) we obtain that M n = 0. Then, the gap equation (2.27) becomes: 



lim 

e^0+ e 



N f N f°° c\k 1 

f = --(N-l) / G n (K) = --(N-l) / —- . (2.30) 

r Jk 9 2 Jo 7T k [exp (k/T) - 1] 

There is no solution to this equation since the integral in the right-hand side is divergent, 
whereas (p 2 — N/g 2 is finite. Thus, the case (ft ^ leads to inconsistencies and cannot 
hold. This means that in two dimensions there is no spontaneous symmetry breaking of 
the global 0(N) symmetry of the nonlinear O(N) model, in agreement with the Mermin- 
Wagner-Coleman theorem mentioned in the Introduction. Note that in four dimensions 
this would have not been the condensation of (ft would take place signalling an 

explicit breaking of chiral symmetry, see for instance Refs. [24-26]. 

(ii) We set (ft = (no spontaneous symmetry breaking). In this case Eq. (2.25) is 
trivially fulfilled. Eq. (2.26) implies that the masses of a and it become degenerate, 

Ml = Ml = M 2 , (2.31) 

so that 

G V {K) = G a (K) = G{K) . (2.32) 

Then, from Eq. (2.27) we find 

g 2 Jk Jo * »l(M) 

where 



u k (M) = \/k 2 + M 2 (2.34) 
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and 

n[u k (M)] = {e W [uj k (M)/T}-l}- 1 . (2.35) 

From the previous equation one can determine the behaviour of the function M(T), see 
later. The case eft = does not lead to any inconsistencies and is the realised one. It 
should be stressed that this solution is highly nontrivial: instead of two massless degrees 
of freedom one has in the vacuum three equally massive particles. 

The basic quantity for a thermodynamic study of the system is the pressure, which is, 
up to a sign, identical to the minimum of the effective potential, 



p = -V e 



mm 
eflf 



From the pressure we can compute the energy density e and the trace anomaly 9 using the 
first principle of thermodynamics: 

e = T dT- p ; 
6 = e — p . 

2.4 Regularisation 

The nonlinear O(N) model is infrared-finite, since the mass is nonvanishing. Thus, there 
are only ultraviolet divergences, which we regularise via an ultraviolet cutoff A [2, 6, 7]: 

where we subtracted divergent M- and T-independent terms. In the previous equations 
the quantity M = M(T) denotes the temperature-dependent mass. At T = one has 
M(T = 0) = m, with m being the only dimensionful parameter of the model (emerging 
through dimensional transmutation) — see the next Subsection for further details. 

Using equations (2.16) and (2.37) we obtain the regularised effective potential and the 
regularised pressure: 

^ NM 2 f A dk k 2 n [u k (M)l NM 2 ( , A 2 \ . 

where the regularised mass M = M(T) is determined from 

N =N r*MnH(M) ] + N ln A? 

g 2 J vr u k (M) 4vr M 2 v ' 



Note that at T = we obtain the following relation: 



2.5 Renormalisation 



After regularisation we can perform the renormalisation. First, we compute the renor- 
malised coupling constant as it was defined in [2]: 
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'o r ~ 

i>ren 



2 d Pieg (T = 0) 
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2 d 
N dm 2 



dm 2 



N- 



. m 



N- 



8^ 



A 2 \ 

1 + In 



1 , A 2 
- — In 

Air fi 2 



(2.41) 



Inserting this result in Eq. (2.38) and subtracting the zero-temperature contribution —Nm 2 / 
we get the renormalised pressure in the limit A — > oo: 



Pi 



NM 2 
NM 2 



25 



ren 



+ A^ 



dfc k 2 n [oj k (M)] NM 2 



7T U) k (M) 8tt 
dk k 2 n [uj k (M )] AM 2 
8~7~ 



7T 



uj k (M) 



1 + ln 



1 + ln 



A 2 \ Nm 2 



M 2 J 
M 2 
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+ 



8tt 
A^m 2 
8vr 



(2.42) 



where we have used Eq. (2.40). The temperature-dependent mass M is determined from 
the renormalised gap equation: 



h- > 
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9ren 


1 


POO 


9ven 


Jo 



1 A 2 
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dfc n [u k (M)} _L_ A2_ 
~ uo k (M) + 47 ""m 2 ' 



dfc n [cj fc (M)\ 1 fj, 2 
it w fc (M) 4vr n M 2 



(2.43) 



Solving Eq. (2.43) at T = one can show that the model exhibits dimensional transmu- 
tation, meaning that there is a nonvanishing mass parameter m in the vacuum which is 
generated due to the renormalisation of quantum corrections: 



■57T 



M 2 (T = 0) 



m 2 = fj 2 exp 



4-7T 
9ren 



(2.44) 
0. Thus, the 



Notice that m depends non-analytically on the coupling constant for g 
mass would vanish if naive perturbation theory were applied. 

In order to determine the asymptotic behaviour of the theory we compute the beta- 
function of the running coupling: 



/'<- 



dg, 



ren 



d/i 
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dyu 
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"2^ 



9' 
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4tt 



In— ^ 

V 2 



4-7T 



ffren 

" 2vr 



<0 . 



(2.45) 



We have shown that the coupling constant decreases with increasing /i and thus the theory 
is asymptotically free. 
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2.6 Two-loop approximation 

The Lagrangian (2.12) contains only three-point interaction terms and no four-point ver- 
tices. Therefore to two- loop order we can only construct sunset diagrams but no double- 
bubbles. Besides, since there is no spontaneous symmetry breaking in two dimensions, 
i.e. cp = and M 2 = M 2 = M 2 , only the term ia(a 2 + 7rf)/2 contributes to the 2PI part 
of the effective potential via a sunset diagram: 



f [ G a (K + P)G(K)G(P) . (2.46) 
Jk Jp 



Taking the sunset diagram into account the effective potential to two-loop order is: 

V cS (a ,G,G a ) = iV^o + ^o (247) 
2g z 2 



+ \ E f K ^Gr\K) + D-r\K;a Q )Gi{K)-l] +V 2 . 



l=cr 1 n,ct 

The corresponding equations for the condensate of the auxiliary field and for the full 
propagators are: 

io » = h{-7 + L G(K) ) (2 - 48) 

G^iK) = D~ l (K) + T,{K) ; (2.49) 

G~\K) = D-\K) + S a (/f) , (2.50) 

where the self-energies for a and for the N scalar fields read 

™ = hmwr L a(p)G * lK ~ p) ' <2 - 51) 

s « (A ' ) = 2 «Sf) = f/ f , G(p)G(A '- p ) • (2 - 52) 

The temperature-dependent mass of the scalar fields is determined by solving 

M 2 = m 2 + ReS = ia + Re J G(P)G a {K - P) , (2.53) 

where we have used m 2 = iato, see Eq. (2.14). Besides, using Eq. (2.15) we obtain the 
following relation for the full propagator of the auxiliary field: 

G- 1 (K)=e + j j G(P)G{K-P) . (2.54) 

Next, we perform the following three steps to simplify these results: 
(i) Out of Eqs. (2.48), (2.53), and (2.54) we obtain 



4+ f G(K) =2e M 2 — Re / G(P)G a (K - P) 
r Jk L Jp 



2e 



M 2 - Re 



G(P) 



J L G(L)G(K-P-L) 



(2.55) 
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In the nonlinear limit e — > + this expression reduces to 
' f G(K) 1 ' 



9 2 Jk J K K 2 + ia + f p G(P)G a (K-P) 

1 

K K 2 + M 2 + Im f p G(P)G a (K - P) ' (2 ' 56) 

where the real part of the sunset has been reabsorbed in the definition of the mass, see 
Eq. (2.53). Moreover, since all three particles have the same mass and decays are strongly 
suppressed, we neglect the imaginary part of the self-energy: 



Im 



j p G(P)G a (K-P) 



. (2.57) 



Then, we obtain at two-loop order the same gap equation as in the one-loop approximation: 

4 = / G(K) ; (2.58) 
T Jk 

this result is consistent with the observation that the next-to-leading-order effective poten- 
tial is minimised by the leading-order mass [6, 7, 10]. 

(ii) Furthermore, using Eqs. (2.15) and (2.54) one can show that in the nonlinear 
limit the two integrals f K D~ l G a (K) - 1 and f K f p G a (K + P)G(K)G(P), involving the 
auxiliary field appearing in the effective potential (2.47), reduce to constant terms which 
are irrelevant for the thermodynamics: 



f D^G^K)-! 

JK 



e + fJ pG (P)G(K-P) e^o + 
Similarly, we get (at the minimum): 

f [ G a (K)G(P)G(K + P) 
Jk Jp 

J G{P)G{K + P) 



N 

T 

N f I 



4 Jk £ + y II G{L)G{K — L) 

N f 1 f 

J G(P)G(K + P) 



£ ->o+ 4 J K § J L G{L)G(K - L 
N f 1 



J G{P)G{K + P) 



l=-q 4 J K f j_ Q G(Q)G(K + Q) 

^ T/ f fjQOW )W + o) X 0(mK + P) ^ onst - (2 59) 

(iii) Finally, we can rewrite the term J K [Z) -1 (.ftr; oto)G(K) — l] as follows: 

' K 2 + m? 



[ [D- 1 (K;a )G(K)-l] = [ 
Jk Jk 



K 2 + M 2 



(2.60) 



(m 2 -M 2 ) [ G(K) 

JK 



K K 2 + M 2 ' Jk X ' 9 



,2 
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where we used the gap equation (2.58). Besides, we set D 1 (K; cto) = K 2 + iao = K 2 + 771 . 

Using the results obtained in (i), (ii) and (iii) the effective potential to two- loop order 
simplifies as follows: 

NM 2 



t /-mm 



NM 2 N f , 1 f -, 

^ + y J^G-\K) + - J K \nG-\K) . (2.61) 

At T = the effective potential can be regularised analytically. Introducing an ultra- 
violet momentum cutoff and dropping m-independent divergences one obtains [2, 6, 7]: 



V eS (T = 0) 



Nm 2 

1 



+ 



+ 



8tt 



Nm 2 ( , A 2 
"5— 1 + ln — 



A 2 A 2 
(A 2 + 2m 2 ) In In - m 2 li 



(2.62) 



+ 



771 

Ait 



7E - 1 - In 



A 2 
4m 2 



with li(x) denoting the logarithmic integral function. 

However, at finite temperature numerical methods must be employed in order to cal- 
culate thermodynamical quantities. As was computed in [6, 7], the final result for the 
renormalised pressure to two-loop order at nonzero T is given by: 



„(2-loop) 



N 



8tt 



M 2 )+N 



dkk 2 n [oj k (Afy)] 

7T U k (M^) 



(2.63) 



NM 2 
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dk n[uj k {M^)] 1 
7T w fc (Afy) 2 



d(F 1 + F 2 ) 
M * dM 2 ~ F '- F2 



Here F\ is calculated using the Abel-Plana formula 



/+00 

and Fi takes the form 

F 2 = V 



dk 
2^ 



lnG-\K 



dk T 



dk 



(27r)7-oo (2tt) 



(2.64) 
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In (k 2 IM 2 



M 2 



■ In In 
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M 2 



where V is indicates the principal value of the integral, and 

1 G-\K) 



4vr y/K 2 (K 2 + 4M 2 ) 



M 2 = M 2 exp 



dk 



/o ^ (M^) exp[w fc (M^)/r]-l 
Besides, AL, is obtained by solving the equations 



(2.65) 

(2.66) 
(2.67) 
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dfc 



1 



1 



A 7 vr w fc (M ) exp [u k (M^/T] - 1 
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+ A 



Ml = M 2 exp 



ln4 + 2vr 



d (Fx + F 2 ) 



dA/ 2 



-I r** 



1 



w fc exp [w fc /T] - 1 



(2.68) 
(2.69) 



- 11 - 



These expressions (which coincide with those of Ref. [6], where a different approximation 
scheme was used) will be used in Sec. 4 when comparing to the lattice results. 



3 Lattice simulation 

In this Section we present the detailed evaluation of thermodynamical quantities of the 0(3) 
model by numerical lattice calculation. The Euclidean, discretised action under study takes 
the form of a Heisenberg model, i.e. an 0(3) spin model with nearest-neighbour interaction: 



where the sum runs over all bonds of a 2-dimensional lattice and the are 3-dimensional 
unit vectors in internal space. The system at finite temperature T is realised by making 
the timelike extent of the lattice finite and consisting of Nt sites, with periodic boundary 
conditions in that direction: denoting with a the lattice spacing, we have aN = 1 /T. The 
identification of the above action with the continuum one defining the partition function 
Eq. (2.1) is completed by setting j3 = N/g 2 . 

In a practical Monte Carlo simulation, moreover, the space extent of the lattice, N x , is 
also necessarily finite; as long as N x 3> Nt, however, this will hardly affect the temperature 
setting. In order to minimise the finite-size corrections, we will assume periodic boundary 
conditions in the spacelike direction throughout all of the following. 

Our goal is to extract the thermodynamical behaviour of the action Eq. (3.1). In 
the next parts, the lattice methods will be illustrated along with the determination of the 
physical scale; subsequently, we present the numerical results that will later be compared 
to the predictions of the previous Section. A discussion of the algorithms implemented, 
their efficiency and autocorrelation is deferred to Appendix A. 

3.1 Thermodynamical observables 

For the evaluation of pressure p, trace anomaly 9 and energy density e, we adapt the 
"integral method" [27] to our (l+l)-dimensional model. We start from the continuum 
expression for the (nonphysical) pressure, 



(the second equality stems from the assumption of spatial isotropy, that always holds in the 
thermodynamic limit), where Z is the partition function associated to S (it is understood 
that it refers to a system with timelike extent Nt): 




(3.1) 



P*(T) = T 



dlnZ 
dV 




(3.2) 




(3.3) 



Since the one-dimensional spatial volume is given by V = aN x , we have 



a 2 p*(T) 



1 



hxZ ; 



(3.4) 



N x N t 
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we then write In Z as an integral in order to express it in terms of measurable quantities, 

a 2 MT)= /"V -i-g^g ! (3.5) 



0o 



N x N t d/3 



0' 



the integrand, exploiting translation invariance, has now become the expectation value of 
the energy per site, 

1 dlnZ /„ \ . , 

x+it), (3.6) 



N x N t 3/3 

with the shorthand i e = Si ■ Si+e- The last step is an additive renormalisation, needed 
to make the pressure finite: the standard choice is p(T) = p*(T) — p*(0), which translates 
into the subtraction of the same integrand evaluated in a zero-temperature system (that 
is, Nt = oo): making coupling and time extent explicit in the average values, we have then 







((4 + ^V^-2(^ i0O )d/3'. (3.7) 



In the above we fix the integration to start at /3q = 0: on the practical side, this value can 
be directly simulated with no particular troubles. 1 Note also that T = will be numerically 
approximated with a large enough square system, that is N x = Nf. 

Likewise, for the trace anomaly, starting from the continuum expression 

* = ^ (£) , (3-8) 



dT \T 2 
one derives 

e ~w = N? (sff ) ( (4 + £t) ^ - 2 ■ (3 - 9) 

This does not require integration over the coupling, but the knowledge of the beta-function 
is needed for the prefactor 

9P dp 
dWT = - a d^- (3 - 10) 
The above will be extracted nonperturbatively, measuring the running of the coupling with 
the scale and using a suitable parametrisation of the resulting data. 
Finally, we also consider the energy density 

T 2 dlnZ „ 

e = T^r = e+ "- < 311 » 

3.2 Scale setting 

The running of the scale a(/3) is needed to convert the temperature in physical units, as 
well as for the calculation of 9 and e. We evaluate it by computing, in a zero-temperature 
system and over a wide range of couplings, the spin-spin correlation function 

C(an) = ■ s i+n e) , (3.12) 



1 It is in fact trivial and corresponds to generating all s; uniformly on the unit 2-sphere. 
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which then, to obtain the zero-temperature mass m, is fitted to the expected functional 
form for a two-dimensional scalar field, valid for large R: 

( e -mR e -m(aN x -R) \ 

C{R) ~ co -= + )+c 1 , (3.13) 

V \JmR yJm(aN x — R) I 

where we consider also the "mirror term" to account for the periodic boundaries. Distances 
on the lattice are in units of a, therefore the masses we measure are the adimensional 
quantities am; once we fix the corresponding physical value we have the function a{(3) 
in units of 1/m. The functional form Eq. (3.13), furthermore, is valid for large enough 
distances: in practice, we will fit the measurements to the function in a variable range 
Rmin < R < Rm&x, looking for a plateau of the resulting am in -R m ; n . Furthermore, as 
the system never undergoes symmetry breaking, we expect the constant c\ to vanish at all 
couplings in the infinite- volume limit. 

3.3 Numerical results 

3.3.1 Beta-function 

The first task is the determination of the running of the scale. We measure the correla- 
tor Eq. (3.12) on lattices of size 156 2 and 256 2 for j3 = 0.1,0.2,- •• ,3.0, generating 10 5 
configurations for each setting, and use Eq. (3.13) to extract am(f3). 

With the system volumes at our disposal, we can fit the data to Eq. (3.13) and find a 
satisfactory plateau for the mass in the region j3 < 1.8: beyond that value, the correlation 
length £ = l/(am) grows too much, making our system size insufficient (Fig. 1, left); 
correspondingly, beyond that value the constant c\ in Eq. (3.13) ceases to be compatible 
with zero. On the other hand, for couplings below f3 ~ 1.0, the correlator signal drops too 
fast, i.e. the lattice becomes too coarse to measure the lightest mass of the theory. 

The evaluation of the coefficient Eq. (3.10) is performed by locally interpolating the 
measured data as 

P(a)=g Q - gi a v (3.14) 

and then differentiating the resulting curve; the interpolation is done on the 256 2 mea- 
surements, since a comparison with the 156 2 data shows that finite-size corrections are 
negligible. In this way we obtain the beta- function shown in Fig. 1 (right). The increasing 
difficulty in fixing the scale as (3 grows will be the main source of systematic error in the 
following results. 

3.3.2 Thermodynamics 

The integrand in Eq. (3.7) is evaluated on lattices with spatial sizes N x = 156 and 256, 
and all time extents Nt = 6, • • ■ ,18 for the finite-temperature part, plus the Nt = N x zero- 
temperature systems. We measure the average energy per site at all couplings ft' < j3 
in steps of 0.1 and perform the integration evaluating the area under the curve 2 , obtain- 
ing the pressure for f3 = 1.0, 1.1, ••• ,1.7, corresponding to lattice spacings in the range 
0.02 <am< 0.55. 

2 The integral was done using a simple trapezoid method. This gives satisfactory results since for the ft 
stepsize which was used the curvature of £(/3) is very small. 
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Figure 1. Left: the running coupling, computed on 156 2 and 256 2 lattices, with two parametri- 
sations to the latter data, Eq. (3.14): fit 1 refers to 0.001 < am < 0.08, while fit 2 is done in 
0.01 < am < 0.6. Error bars are much smaller than the symbol size; only some sample points are 
shown for the 156 2 lattice. Right: numerical results for the beta- function. Full circles indicate the 
final values used to set the physical scale. 
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Figure 2. The physical pressure p(T/m)/T 2 from N x — 256 systems at various lattice spacings. 
Each symbol represents a different (/3, N x ). Error bars show the systematic error from both discreti- 
sation and finite-size effects. The lines represent analytical calculations at LO and NLO (Eqs. (2.42) 
and (2.64), respectively); the horizontal line marks the asymptotic value ^. 



For each combination of Nf,N x and (3, 10 5 measurements were taken. The statistical 
error is estimated by binning the data in chunks of 1000, using Eq. (3.7) and averaging the 
resulting points for the physical pressure: it turns out that the uncertainties are very small 
and completely overcome by systematics. 



- 15 - 



0.4 
0.35 

0.3 
0.25 

eg 

t; 0.2 

<x> 

0.15 
0.1 
0.05 


1 2 3 4 5 6 7 

T/m 

Figure 3. Trace anomaly 9(T/m)/T 2 ; symbols are as in Fig. 2. 

The results for p are plotted in Fig. 2, along with the analytical expectations. The slight 
differences between discretisations and system volumes are used to assess the systematic 
errors; however, the neat agreement between the various data sets is such that we consider 
both the thermodynamic and continuum limits to be effectively reached. 

Similar considerations hold for the trace anomaly, calculated using the same input data 
through Eq. (3.9), and the energy density, obtained with Eq. (3.11). These are plotted in 
Figs. 3 and 4 respectively, along with the analytical predictions. 

4 Comparison of analytical and lattice results 

In this Section we compare the analytical results obtained in Sec. 2 with the lattice ones 
computed in Sec. 3 for the case N = 3. We organise this Section by first discussing 
pressure, trace anomaly, and energy density; then, we turn to two consequences of our 
study: the behaviour of the quasi-particle mass at nonzero temperature and the role of 
instantons/calorons in the thermodynamics. 

Pressure. In Fig. 2 the pressure is plotted as function of T: we present the analytic 
curves at one- and two-loop order and the lattice points for different values of j3. As one 
can deduce from the plot, the one-loop approximation describes well the lattice data for 
small temperatures. However, it overshoots them for large T. The pressure at two-loop 
level considerably improves the agreement of theory and lattice simulations for high T. 
Although a departure is still visible for increasing temperature, we can conclude that there 
is a good and parameter-free agreement between these two very different methods (the 
CJT formalism and lattice) of evaluating the pressure of the system. 




- 16 - 



1.2 



1 - 




1 2 3 4 5 6 7 

T/m 



Figure 4. Energy density e(T/m)/T 2 ; symbols are as in Fig. 2. 

Notice that this finding could be expected by a closer inspection of the properties of 
our analytic expressions: in fact, at one- loop level a good agreement with lattice data is 
not achievable because of a wrong asymptotic limit. Namely, in both cases the pressure 
approaches for high T the limit of a non-interacting gas of bosons, where each degree 
of freedom contributes as ttT 2 /6. However, at one-loop level there are, at high T, N 
active degrees of freedom, thus p( 1_lo °P)(T 3> m) = NirT 2 /6 for T — > oo. This cannot 
be correct because the nonlinear constraint in Eq. (2.4) eliminates one degree of freedom. 
On the contrary, at the two-loop level one obtains the correct result p( 2_lo °p)(T 3> m) = 
(N — 1)ttT 2 /6 for T — > oo, which shows that only N — 1 non-interacting bosons are present, 
as it should be. This fact explains the much better agreement of the two-loop approximation 
at high temperature. 

Trace anomaly. In Fig. 3 we present a similar plot of the function 6/T 2 . The apparent 
larger mismatch at high T relies on the fact that 6 = e — p = T 3 8t{p/T 2 ): the presence 
of a derivative enhances the disagreement for large T. However, it is very interesting to 
stress that both approaches show a form of the function 6/T 2 which looks very much like 
the one obtained in the deconfined phase of Yang-Mills theory, e.g. [27, 28]. In fact, the 
present two-dimensional model, which does not have any phase transition, represents a 
simplified description of a gluonic gas above the critical temperature. Note that also in 
lattice simulations of Yang-Mills theories the maximum of the trace anomaly is reached 
above T c just as our low-dimensional model does. 

Energy density. In Fig. 4 we present the plot for the energy density. Analogous 
comments regarding the comparison of lattice and analytic results hold. Moreover, here it 
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Figure 5. The quasi-particle mass M(T) as function of T. For T — > the constant value M(0) = m 
is reached (dimensional transmutation). For high T the perturbative behaviour M{T) ~ T/ln(T) 
is realised. 



is clearly visible that the one-loop result overshoots the corresponding correct Boltzmann 
limit of (N — 1) active DOFs in the high-T regime. 

Quasi-particle mass. In the framework of our analytic approach we have calculated 
the mass of the particles as function of T. In Fig. 5 we plot the mass of the scalar particle 
M(T) normalised to its vacuum value M(0) = m as function of the temperature (notice 
that the mass M(T) is the same for both the one-loop and the two-loop calculations). 
The function M(T) starts from a nonzero value at T = and then bends and reaches, as 
expected by basic dimensional consideration, a straight line for large T. This behaviour 
is reminiscent of that of the gluon mass in the deconfined phase of Yang-Mills theories 
[29]. The model exhibits dimensional transmutation just as QCD, meaning that at zero 
temperature there is a nonvanishing mass gap, which is generated due to renormalisation 
of quantum corrections. Note that the temperature at which the function M(T) bends 
corresponds roughly to the maximum of the trace anomaly. Conversely, at high T the 
temperature-dependence of the mass can be approximately parametrised by T/logT in 
accordance with the aforementioned results concerning Yang-Mills theories. 

Role of instantons. The 0(3) model in 1 + 1 dimensions exhibits instantons in the 
vacuum and calorons at nonzero temperature. It is natural to ask if and how these topo- 
logical objects affect thermodynamics. To this end, it should be stressed that the analytical 
method presented in Sec. 2 does not distinguish between trivial and nontrivial filed config- 
urations, neither in the vacuum nor at nonzero temperature. At a first sight, one may then 
argue that no instanton effects are present in our model: in fact, the method employed is 
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valid for each N, while instantons are a peculiarity of the N = 3 case only. However, a 
closer inspection of our results is needed: namely, the system — as illustrated in Sec. 2 - 
shows nonperturbative properties: the vacuum is not realised by a constant value of the 
fields, such as a = <p ^ 0, rather by <j) = 0. Considering that the simple field configura- 
tions g = 0, 7Tj = is obviously not allowed (it violates the constraint), this means that 
rather complicated functions for a(x,r) and 7Tj(x,r) with zero average are realised. For 
instance, an ensemble of calorons and anticalorons fulfills the required properties: all fields 
have zero average but nonzero contributions. We thus argue that the role of topological 
objects is prominent, because it naturally explains the nontrivial emergence of the vacuum 
(a) = <7Ti) = 0. 

Indeed, in Ref. [21] the 0(3) model has been compared to the equivalent CP 1 model. 
In the case of CP 11-1 models a calculation for the pressure was done only in the topologi- 
cally trivial sector: in fact, in that case it was possible to distinguish between the different 
topological sectors of the theory. A comparison of the pressure of the CP 1 model with 
that of the 0(3) model shows a large discrepancy, which was ascribed to instanton/caloron 
effects by the authors of Ref. [21]. Obviously, this conclusion holds only under the as- 
sumption that the instanton contributions are indeed included in the thermodynamical 
treatment of the 0(3) model. As we argue above, this seems to be the case. However, 
a direct calculation of the instanton role in the 0(3) model would be very desirable from 
both the analytic and the lattice sides. Moreover, the 0(3) calorons contain also monopole 
constituents, just as Yang-Mills theories [15]. An estimate of the contributions of these 
constituents to thermodynamics would also be very interesting. 

5 Summary 

In this work we have studied the thermodynamics of the 0(N = 3) model in 1+1 dimensions 
using the CJT formalism and the auxiliary field method to one- and two-loop orders. At 
the same time, we have performed a precise lattice simulation for this system at nonzero 
temperature applying the integral method. 

Our main results are displayed in Figs. 2, 3 and 4 where the theoretical curves have 
been compared to the lattice calculation. It should be stressed that, once the mass of 
the field in the vacuum is used as an energy unit, no free parameters in the analytical 
calculation are present: the agreement shown in the aforementioned figures confirms the 
power of the employed analytic formalism. In particular, at one-loop level the results are 
very good for small temperature, but not so at large T because of an expected failure 
in this energy domain. Namely, in this approximation a gas of free N bosons is realised 
in this limit, while only N — 1 degrees of freedom exist due to the nonlinear constraint 
entering the definition of the model. At two-loop order this problem is solved and the 
correct asymptotic limit is reached. The agreement of theory and lattice is therefore much 
better in the high-T domain. Surely, theoretical improvement can still be done in order to 
achieve an even better matching. 

The O(N) model is very interesting because it shares some properties with Yang-Mills 
theory, such as trace anomaly and asymptotic freedom. In fact, the form of the trace 
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anomaly (see Fig. 3) shows a striking analogy to the more difficult four-dimensional Yang- 
Mills case. For these reasons the lattice results presented here can be of interest for testing 
other theoretical methods dealing with nonperturbative features of quantum field theory. 
Namely, while going beyond the two-loop CJT results presented in this work would be a 
very hard technical task, one could develop new theoretical methods which involve new 
resummation schemes. Surely, a better understanding of such methods in a model which 
shares many similarities with Yang- Mills gauge systems can be of great use to improve our 
knowledge on the latter. 

A further interesting achievement has been the determination of the temperature- 
dependent mass M(T), see Fig. 5. Also in this case a form which is reminiscent of that of 
the gluon mass in the deconfined phase of Yang- Mills systems is obtained: for small T one 
has M(T) ~ m 7^ (dimensional transmutation), while for large T one has the expected 
perturbative growth M(T) ~ T/lnT. 

For the particular choice N = 3 used in this work the 0(3) model contains also 
instantons. We have argued that these topological configurations, although not explic- 
itly included in our evaluation, are anyhow present because they naturally generate the 
nontrivial vacuum of the system. Indeed, a further study in this direction represents an 
interesting outlook both as an analytic and as a lattice study of the 0(3) system. 

Finally, one can also include in the present model a so-called #-term which mimics 
the axial anomaly of QCD (see for instance [30]): a thermodynamical investigation of this 
system is feasible both analytically and numerically. 
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A Algorithms and autocorrelations 

We now come to a study and comparison of the various algorithms that can be imple- 
mented to investigate numerically the model, Eq. (3.1). Aside from the ordinary, quite 
inefficient Metropolis algorithm, which we use only for cross-checking purposes, we tested 
the heatbath update method, the cluster method and the microcanonical overrelaxation. 

In the following we measure and compare the most important indicator of an algo- 
rithm's efficiency, namely the (integrated) autocorrelation time Ti n t, which measures how 
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many elementary update steps are needed to generate a new, statistically uncorrelated con- 
figuration while exploring the phase space of the theory through the Markov chain driven 
by the update method. Furthermore, the final measure of the efficiency of an algorithm is 
not Ti n t itself (measured in numbers of updates), but rather -rf^ 1 (measured in seconds), 
the product of T{ n t and the average wallclock time it takes, on a given machine, to perform 
a single update (which, in particular for cluster-based methods, may largely depend on the 
coupling). All simulations presented in this work have been done on individual cores of 
AMD Opteron 2427 (Istanbul) CPUs with a clock rate of 2211 MHz. No parallelisation 
has been employed. 

In general, different autocorrelations are associated to different observables. We focus 
on three quantities: two are of obvious physical interest, corresponding to the energy- and 
ordering-like relevant scaling operators, which are general for O(N) models, 

E< W )(i-3-3f)) 

m = x ' . (A.l) 



2N x N t ' N x N t 

The third one is the volume-average of a single component of the Sj and gives a measure 
for "isotropy-like" correlations, i.e. a memory of prior global direction in internal space, 



(A.2) 



N x N t 

Its autocorrelation, which is technically interesting and can in fact become quite large 
(where the system has large spatial correlations), is however for the most part irrele- 
vant from the quantum field theory perspective, which is concerned primarily with 0(3)- 
symmetric physics. 

To estimate the autocorrelation, one can employ a binning analysis with bins of variable 
size [31], or use the jackknife method; here, we apply the more refined technique proposed 
in [32]. All autocorrelation measurements have to be made after thermalisation is reached, 
which in the present case is determined by inspection of the Monte Carlo time histories for 
the observables defined above. 

In the following we discuss the algorithms we implemented in terms of their Tint and 
7j^, and present the mixed strategy we adopted for all simulations reported in Subsection 
3.3, which ensures efficiency over the whole /3-range of interest. 

Heatbath 

The model, Eq. (3.1), is suited for the heatbath update. The new vector on site i is 
extracted randomly according to the equilibrium probability distribution induced by the 
associated contribution to the action, involving its four neighbours: 



P(si) ~ exp fisi • ^2 Si+v (A.3) 



In practice, it amounts to generating the angle between the sum vector and the new Si 
according to the above formula, and then uniformly choosing a random azimuthal angle 
[31]. 
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Figure 6. Real-time t[^ 1 from simulations on a 64 2 lattice, for heatbath (left), Wolff (middle) and 
mixed (right) update methods. Note the semi-logarithmic scale. 

All sites of the lattice are updated in a checkerboard fashion: this constitutes an 
update step ("sweep"). In Fig. 6 (left) the behaviour of t^ 1 for the three observables 
on a 64 2 lattice is shown (we focus on this quantity instead of T m t to make a meaningful 
comparison between different update methods, however for the heatbath update the two 
are proportional in a /3-independent way). It is apparent that the method, while performing 
well for small couplings, suffers from the long-range correlations that develop as > 1, 
especially for the ordering-like operators: there, alternative methods should be considered. 

Cluster algorithm 

We then test the cluster update method, developed by Wolff [34] generalising the work of 
Swendsen and Wang [33] for the Ising model. The idea is to build clusters of spins, reflecting 
to some extent the structures, in the system, which carry the long-range correlations, and 
update these clusters all at once. For our model, an update step consists of the following 
operations: 

a) Choose a random site Xi and a random reflection plane tt with surface normal n 
("seeding"). Reflect s~l across tt and mark 2j. 

b) Examine all nearest neighbours {xj} of xi and with probability 

p = 1 — exp{min[0, 2/3(n • sl)(n ■ s~j)]} 
reflect s5 across ir and mark well. 

c) Repeat step b) for all links connecting newly marked sites to unmarked neighbours 
until the process stops. 

Clustering methods are generally understood to be free of critical slowing down: thus, 
they are suitable for simulations in the vicinity of a critical point. Our model does not 
exhibit a true transition in two dimensions, nevertheless we implement the Wolff algorithm 
and investigate its efficiency, since on any finite volume we see a strong ordering tendency 
at large enough couplings. 

The average size of a cluster (therefore, the time needed perform a single update) grows 
with f3: indeed, we observe it increases by more than an order of magnitude between f3 = 0.1 
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Figure 7. Tint as a function of the coupling for the mixed update on square systems of various 
sizes, for the e (left) and the m (right) observables. Note the semi-logarithmic scale. 

and (3 = 3.0. This is, however, completely countered by a shortening of the autocorrelation 
as the clusters span a larger portion of the system, so that, in terms of T^ t al , the efficiency 
of this algorithm in disordering the configuration is generally better at larger couplings 
(Fig. 6, middle). In the region of small (3, the heatbath is still superior than this method, 
so in the following we devise a mixed update that ensures small autocorrelation over the 
whole range of couplings under study. 

The mixed strategy 

Based on the findings presented so far, we expect that a judicious combination of the two 
algorithms would perform well at all couplings. We then define from now on an "update 
sweep" as the combination of a single heatbath sweep throughout the lattice combined 
with ten cluster seedings. 

Fig. 6 (right) shows the resulting r^ t al on a 64 2 system: this update is indeed more 
efficient than either heatbath or cluster alone over the whole j3 range. Hence, we use the 
mixed strategy for the thermodynamical investigation presented in Subsection 3.3. 

We now turn to the question of the scaling of r; nt with the system volume for the mixed 
update: this is useful for a proper tuning of the Monte Carlo data taking based on the size 
of the lattice under study. We measure Tint m units of sweeps, as defined above, on various 
square systems with size up to N x = 64 for energy- and ordering-like observables; these 
are plotted in Fig. 7. The latter seems always to be a bit slower to decorrelate, so we base 
our tuning on its behaviour. It turns out that the height peak of T{ n t(f3), located roughly 
at f3 ~ 1.2, can be described as a function of the system size with a power law: 

Tha*(N x ) ~ aNf , a = 0.13(5) , b = 0.38(6) , (A.4) 

this implies the upper bounds Tint* (156) ~ 6 and Tj nt *(256) ~ 9; to be on the safe side, in 
collecting the data of Subsection 3.3, we perform respectively 14 and 21 sweeps between 
measurements. 

Metropolis and overrelaxation 
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For the purpose of cross-checking, we also implemented a standard Metropolis update, 
where a proposed new value of s*j is accepted with probability p = min[l, exp(— AS 1 )], with 
AS the corresponding change in the action. The associated autocorrelation is qualitatively 
similar to the heatbath procedure, but it is about an order of magniture larger than the 
latter. 

We also investigated whether the algorithm's efficiency can be improved with micro- 
canonical overrelaxation, i.e. a random rotation of about the sum of the neighbours' 
spins so to achieve maximum disordering of the configuration while exactly preserving the 
energy and the Metropolis acceptance. While this brings indeed the Metropolis update to 
an efficiency comparable to the heatbath in the region /3 > 1, we do not see any gain at 
lower couplings: therefore, we refrain from exploring this avenue any further. 
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